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ABSTRACT 



Apparatus for transmitting and receiving digital information 
over mobile communication channels. The apparatus 
includes a pilot symbol insertion device for periodically 
inserting data dependent pilot symbols into frames of digital 
data for subsequent channel estimation, a continuous phase 
modulator for modulating and transmitting over a mobile 
communication channel the frames of digital data and means 
for receiving and filtering the transmitted data. A channel 
estimator estimates the channel amplitude and phase distor- 
tion from the received pilot signal at different time instants. 
Means are provided for demodulating the received signal 
and for decoding the demodulated signal to recover the 
estimated transmitted data. 

29 Claims, 9 Drawing Sheets 
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PILOT-SYMBOL AIDED CONTINUOUS 
PHASE MODULATION SYSTEM 

FIELD 

The present invention relates to a method and apparatus 
which uses data dependent pilot symbols in continuous 
phase modulation (CPM) systems for the transmission and 
detection of digital data in mobile radio applications to 
reduce the bit error rate of the CPM. 

BACKGROUND 

There are different types of interference that can lead to 
reduced detected signal strength of a transmitted signal at a 
receiver. Reflections or multiple sources which result in 
phase differences at the receiver can reduce the detected 
signal strength. Another type of interference which reduces 
detected signal strength is known as flat fading, in which the 
frequency response of the received signal is flat and only the 
gain and phase fluctuate. 

The use of pilot symbol assisted modulation is a technique 
that has been applied for nunimizing the effect of flat fading, 
and is particularly useful with mobile receivers. However, 
the use of pilot symbols for channel estimation has been 
confined to linear modulation schemes. It has been generally 
regarded that accurately estimating the phase of the received 
signal in a coherent detection of CPM signals is difficult and 
so little attention has been given to this problem The 
difficulty arises because the transmitted CPM signal in any 
given interval depends not only on the current input symbols 
as well as the past input symbols. Thus, conventional pilot 
symbol techniques which are data independent and known 
channel estimation techniques can not be applied to CPM. 

SUMMARY OF THE INVENTION 

According to the invention there is provided an apparatus 
for transmitting and receiving digital information over 
mobile communication channels. The apparatus includes a 
pilot symbol insertion device for periodically inserting data 
dependent pilot symbols into frames of digital data for 
subsequent channel estimation, a continuous phase modu- 
lator for modulating and transmitting over a mobile com- 
munication channel the frames of digital data and means for 
receiving and filtering the transmitted data. A channel esti- 
mator estimates the channel amplitude and phase distortion 
from the received pilot signal at different time instants. 
Means are provided for demodulating the received signal 
and for decoding the demodulated signal to recover the 
estimated transmitted data. 

Advantageously, there is provided a pre coder employing 
a finite state preceding method whose operation is based 
upon a trellis structure substantially the same as the con- 
tinuous phase modulation scheme. 

The trellis structure may be that of a full response CPM 
scheme with an identical number of modulation levels and 
modulation index to that of the underlying CPM scheme. 

The coherent detection technique for Continuous Phase 
Modulation (CPM) operating in the Rayleigh flat fading 
channel is based on the idea of inserting periodically data 
dependent pilot symbols that force the CPM signal to pass 
through known phase states. This transmission format 
enables the receiver to extract from the received signal the 
channel fading gains at regularly spaced instants. When 
coupled with proper channel estimation filters, very accurate 
channel state information (CSI) can be estimated at the 
receiver for fading compensation. Moreover, the accuracy of 
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the CSI can be further refined by adopting a multiple-pass 
decoding approach. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 The novel features believed characteristic of the invention 
are set forth in the appended claims. The invention itself, 
however, as well as other features and advantages thereof, 
will be best understood by reference to the description which 
follows read in conjunction with the accompanying draw- 
io ings wherein: 

FIG. 1 is a block diagram of the CPM system; 
FIG. 2 is a block diagram of the encoder; 
FIG. 3 is a trellis diagram of a binary precoder; 
15 FIG. 4 is a trellis diagram of a 4-level precoder; 
FIG. 5 is a block diagram of the channel estimator; 
FIG. 6 is a schematic diagram of a frame showing the 
position of the pilot symbols; 

FIG. 7 is a graph showing the analytical and simulation bit 
20 error rate results for Minimum Shift Keying (MSK) at a 
normalized fade rate of 0.01 with precoding; 

FIG. 8 is a graph showing the analytical and simulation bit 
error rate results for Minimum Shift Keying (MSK) at a 
normalized fade rate of 0.03 with precoding; 
25 FIG. 9 is a graph showing the analytical and simulation bit 
error rate results for Gaussian Minimum Shift Keying 
(GMSK) at a normalized fade rate of 0.01 with precoding; 

FIG. 10 is a graph showing the analytical and simulation 
bit error rate results for Gaussian Minimum Shift Keying 
(GMSK) at a normalized fade rate of 0.03 with precoding; 
and 

FIG. 11 is a graph showing the effect of precoding in 4 
level Continuous Phase Frequency Shift Keying (4 CPFSK). 

35 DETAILED DESCRIPTION WITH REFERENCE 
TO THE DRAWINGS 

The block diagram of the CPM system, shown in FIG. 1. 
includes an encoder 10 which receives the input data which 
40 is a sequence of binary Q-tuples b=( . . . , b_|, b©. bj, — ). 

where b k =Q>\ b 2 b°), Q=log 2 M. b m e{0.1}, and M is 

the number of modulation levels. The encoder divides the 
sequence b into frames of K consecutive Q-tuples. Each 
frame is transformed by the encoder 10 into N channel 
45 symbols; N>K The sequence of channel symbols is denoted 
by c=( . . . . c_ £ , Cq, c_i. . . . ), where each c k is chosen from 

the set {±1, ±3 ±(M-1)} with equal probability. Once 

all the data bits within the same frame have been precoded 
into channel symbol, pilot symbols are then added. 
50 The channel symbols are output to a CPM Modulator 12. 
The CPM Modulator 12 carries out a modulation scheme 
which may have two or more levels. The resultant modu- 
lated signals s(t) are sent over a mobile radio channel with 
a Rayleigh flat fading model. The symbols g(t) and n(t) 
55 represent a zero mean complex Gaussian random process 
and represents flat fading and n(t) represents the channel's 
additive Gaussian noise. The transmitted signal is received 
by an ideal anti-aliasing filter 14 and the detected signal r(t) 
is filtered with a one-sided bandwidth of f J2 Hz. The filtered 
60 signal is sampled by a sampler 15 and fed to a channel 
estimator 16 and a Viterbi Demodulator 18. The channel 
estimator 16 estimates from r the channel fading process 
g(t). This channel state information (CSI) is men exploited 
by a Viterbi receiver 18 during the coherent demodulation of 
65 the CPM signal The output of the Viterbi Demodulator 18 
is sent to a decoder 20 which undoes the mapping done by 
the encoder 10. 
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The structure of the encoder 10 is seen in FIG. 2 consists 
of a precoder 22 and a pilot symbol insertion device 24 

— L-l f tt (3) 

in conventional systems, a memoryless one-to-one map- £o qiiT> = J 1 

per is used as the precoder. However, precoder 22 is a trellis 5 ° 

precoder which reduces (relative to a memoryless mapper) ^ addition, only those f(t) values that are symmetrical about 

the decoded bit error rate in M-level CPM systems with a m included - F <* these pulses, 

rational modulation index of h=p/M. Examples of a binary, ^ JTi 

h=^ and a 4-level, h=V* precoder trellis stmcture are shown ^ <«> 

L"P S ' 4 res P ectivc, y- Referring to FIG, 3 if the » Note that pulse symmetry is very important in arriving at 

inittal precoder state is "0" on the top left and goes to a "1" efficient pilot encoding rules. 8 * 

state it produces a channel symbol of "1" as a result The The term ^ in (1) can be expressed recursively according 

initial precoder state must be even. The next precoder state ^ the foUowing equation: J B 

must be odd. If the next change is from the "1" state back 

to a "<r state then a "-r channel symbol is generated. The 15 fc^^V^T) (5) 

two precoder structures are matched to error event charac- *=° 

tensbcs of the underlying CPM schemes. In each case, the When taken modulo 2ju <f> becomes the phase state of rh. 

precoder 22 only introduces 1 bit of error to the most modulator 12 at time nT L TtaSSL ^L^^lj! 
S^hT" 20 ^**antisobtaine^^ 

introduced by a memoryless mapper. This factor of 2 reduc- ■ • « <W**> with the phase sta*Tu " ^ 

uon in the number of erroneous bits implies a potential 3 dB 

gam in power efficiency in Rayleigh fading channels with an S " H *~ c " c — ' ■ - ■ • (6) 

inverse signal-to-noise ratio error characteristic. Tt ,- . ^ . a . . 

Implicit m the definition in (6) is that for full response CPM 
Once all the data bits within the same frame have been 25 «■* modulator state is the phase state. The coUection of the 

precoded into channel symbols, pilot symbols are then s y stem ***** ™ d *e corresponding state transitions from 

added by pilot symbol insertion device 24. As in pilot one ***** mstant to *e next forms the CPM trellis 

symbol assisted linear modulations, the added symbols are ^ ^^mitted CPM signal, in complex baseband 
required for accurate channel estimation. However because ^ notatlon - * 

of the inherent memory in CPM. the pilot symbols being _ A ~ 0 

inserted are dependent on the input data. The rules for (7) 

generating pilot symbols for M-level CPM schemes with a where A is the transmitted signal amplitude Thi, • 

mc^tionindex of 1/M are given later. The generalization sent over a mobile ramol^ 

LmI^L^ 06 10 h=lVM ' Where pisa » ^ger 35 model. The corresponding recete^ j£R ^ 8 
and M the number of modulation levels, is also given later 

It will be shown that all full response CPM schemes, as well w 

as those partial response schemes using Nyqirist-3 pulse The t«m «*v. ; 

shaping (but without correlative encoding), require oiUvone n^J^r^V * ^ "J™ f 0mplex random 

of channels and K is the number of consecutive Q-tuples per "wcnon ror g(t> is 

frame. For general partial response schemes such as gener- ^^^Wt^o/JoC^x) ^ 

akzed tamed frequency modulation (GTFM) and Gaussian u 2 • 

Minimum Shift Keying (GMSK). more than one pilot sym- ^ T * 15 &C average P° wer * e fading process, f 0 is 

bol is required per frame. 45 ™ e maximum Doppler frequency, J 0 (*) is the zero-th order 

^/T^ -Vj • • c -<- • • • > of — ---^ n ^^fa 

channel symbols arc fed to the modulator 12 at a rate of one idealized scattering model, it is used for the rXose 

symbol perT second in a frame as shown in FIG. 6 where of developing analytical performance results. A properly 

and Q -are : Pilot Symbols. The phase trajectory of the » truncated version of R,(t) will be used for generating 
resulting CPM signal in the interval [(n-l)TjiT]. where n is simulation results. The function n^t) is the channel's addi- 
an integer, is given by t»ve white Gaussian noise (AWGN). It is also a zero mean 

C ^ ,e ^L GaUSsian process and has a power spectral density 

«(0=^, + «*g^-»r + n <'> sj delink (SNR) is 

wh^((,„ismevalueof((<t)attin»enT.hismeiiKKlulation E tr . (i0) 

index, and q^t) is a waveform derived from the frequency = *W d ' = 

pulse f(t) according to the following: 0 J 0 N « N « 

60 Since each frame of QK data bits is encoded into N channel 



^ J o<,<, <2> ^^.^^^^.E^ is defined as 

I o ° * Conceptually, we can assume that the received signal r(t) 

1 01heTOBe ^ fi^t Rsses through an ideal anti-aliasing filter with a one- 

fi is assumed that f(t) is strictly time-limited to the interval „f ^ " 2 ^ mtered ^ " ""'P'ed at 

.0. IT, and ma, its area is normaJized to unity. This means: frequency isTtilst £ud £TfX 
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integer. In the absence of adjacent channel interference CPM scheme with h=l/M is identical to that of 7tfM— 

(Ad), the sampling rate can. in theory, be as large as MPSK. Specifically, there are 2M states, labelled M (T to 

possible. Consequently the intersymbol interference (EI) "2M-F\ The modulator 12 can change state from any of the 

generated from the anti-aliasing filter can be ignored. In the even states to any of the odd states, and vice versa, 
presence of ACL the sampling rate f t must be chosen so that 5 Consequently, there are always M transitions emanating 

the total ISI and ACI is minimized. from and merging into each state. In the transition from state 

The output samples of the anti-aliasing filter is denoted by *y at time (n-l)T to state T at time nT. the channel symbol 

r=( r_ x . r Q , r, ). where c ^ satisfies: 

™* n * (U> 10 A -M<A< M (»> 

g, is the value of the fading process g(t) at time t=CTV s, is Cm = A-2M M<A<2W 

the transmitted signal s(t) evaluated at the same time instant " & + -2M <&<m 
and n, is a sample of the filtered Gaussian noise. The 

different n/s are statistically independent and identically where A=i-j is always an odd number in the range 1-2M. 

distributed (iid), each having a zero mean and a variance of is 2M], Note that transitions emanating from any given state 

a2^f^ a . The sequence r is passed to the channel estimator rece j V e all the symbols from the signal alphabet {±1. ±3. . 

16 shown in more detail in FIG. 3. The latter consists of both ,±(M-1)}. The same is true for the transitions joining at 

a front end filter 25 and an interpolation filter 26. The ^ y . eQ state 

function of the channel estimator 16 is to estimate from r the * ' ^ti; e ctr,.H..» f w 

channel fading process g(t). This channel state information 20 Now if for the precoder the same trel^ s^cn^e (see 

(CSI) is then Exploited by a receiver of the Viterbi demodu- FIGS. 3 and 4) and the same encoding rule (13) is used as 

lator 18 during the coherent demodulation of the CRM for the modulator 12. then the design of the precoder 22 is 

signal The channel estimator 16 and the Viterbi demodu- equivalent to finding the optimal state dependent assign- 

lator 18 are discussed in more detail later. ments of input binary data to the transitions in the modula- 

The Precoder 25 tion trellis that minimizes the decoded bit error rate. In this 

The output of the Viterbi demodulator 16 in FIG. 1 be study, we adopt the heuristic rules that (1) all the transitions 

denoted by c^ Like the transmitted symbol c*. each detected joining at any given state have the same input binary pattern, 

symbol c. takes on a value from the set {±1. ±3. . . . and (2) Gray mapping is used to assign input binary patterns 

4<M-1)}. These symbols are fed to the decoder 20 that to all the transitions diverging from any given state. These 

undoes the mapping performed by the precoder 22. so rules guarantee that the most dominant error events always 

For binary CPM schemes with a modulation index of receive only one bit of error, as oppose to the 2 bit of errors 

h=V*. an error event in the Viterbi demodulator 16 can found in a Gray mapper. However, it should be pointed out 

always be described in the form that the precooVrs 22 in FIG. 2 actually introduce more 

errors to longer error events than the memoryless mappers. 

r=<Yi. • • - M+2A0, . ...0A+2A0,. •) (12) 35 Fortunately, these longer error events occur much less 

where Y*=c*-c* is called a difference symboL and without frequent Consequently, me precoders in FIG. 2 are able to 

loss of generality we assume Yi is positive. Note that for an provide a net performance improvement 

error event of length 1 in the demodulation trellis, there are Finally, it should be pointed out that the preceding rule 

altogether 1-1^ number of "0 M symbols in between the given by (13) is only valid for h=l/M. For h=pfM. where p 

difference symbols "+2" and 4 *±2". where 1^ is the length 40 i s an integer, a different preceding rule can be derived (once 

of the shortest error events. This error event structure again) from the CPM trellis. Although the precoders dis- 

depicted by (12) means that when a one-to-one memoryless cussed are derived from the trellises of full response CRM 

mapper is used to encode the actual data into channel schemes, they can be applied equally well in partial response 

symbols, the binary CRM demodulator always introduces 2 systems. This stems from the fact that both full and partial 

bit of errors to an error event independent of its length. 45 response CPM schemes have identical error event structures. 

However, if the precoder 22 in FIG. 2 is used instead, mere The Pilot Symbol Encoding Rules 

will be one bit of error in the initial split no error in the final ^ shown iously me Mpat of ^ CPM modulator 12 

remerge. and 1 bit of error for each of the ^symbols in ^ depends not only on the current channel 

between the difference symbols +2 and I ±2 JJince the s but ^so on the phase state as well as the previous 

performance of a CPM scheme is dominated by- the shortest 50 ^ Because of this memory effect the pilot 

error events, this means the precoder in FIG. 2b can provide s ^^ x eocoding ^ used in Unear modulations can no 

up to 3 dB of coding gain in fading channels. a ^ applied to CPM. The following are data 

The structure of the error events for M-level CPM dependent pilot encoding rules for M-level CPM schemes 

schemes can be found in equation. 3 of [2J. In the case of ^ fa=1/M Tfaese nlles guarante e that the modulator can be 

M=4 and h=«/4. the shortest and the most^dominant error 55 ^ rcturn riodically to a known phase state. The 

events are of the form T=<+2.-2.0. . ), or IM-M5.--6.0. . . . ), generalization of the proposed pilot symbol encoding tech- 

or TM-N3.+2.0. . . . ). Without loss of generahty.it is assumed * ^ fa an mte ^ fa ^ ven ^ For convenience. 

Yi is positive. With the precoder shown in FIG. 2 each of ^ ^ *^ ^ b chosen ^ the one which me modulator 

these error events causes only one bit of error, which is half Uis forced to return to periodically. For some appUcations 

that introduced by a memoryless Gray inapper. Sirnilarly. for 60 it would ^ desirable to randomize these known states so that 

all the remaining shortest error events, which are of the form e moduUted signal will ^ contain any spectral lines. 

IM+4.±4,0 ). the precoder 22 in FIG. 2 introduces 2 bits vTl «. 

of error while a Gray mapper introduces 4 bits of error. Thus Ful1 Response Signaling 

once again, it can be concluded that precoding can substan- Considering first the case of full response signalling m 

tially reduce the bit error rate of the system. 65 which the pulse length, L. equals unity. Assuming an initial 

The design of the 2 precoders mentioned above is based phase state of <|> 0 =0, then repeated use of (5) enables us to 

on the observation that the trellis diagram of a full response write <)>„ as 
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4w = jiA X o mod 2rr 
If h=l/M and n is odd, then <|> n must be of the form: 
icmod27c;me{±l > ±3, . . . , ztfltf - 1)}, 



(14) 



this section can be generalized to rational modulation indi- 
ces of the form h=p/M. 



5 as 



By setting the term <$> 0 to zero, equation (5) can be written 



(15) 



where the convention that all the phase states are defined 
over the interval of [-it, n] is adopted. In this case, one can 
torce 4> A to zero by using a pilot symbol c^^-m, which 
is a valid channel symbol. On the other hand if h*l/M, it can 
be shown that one cannot, in general, force 6^., to zero in 
one-step by using a valid channel symbol from the set {±1, 
±3, . . . ^M-l)}. If the restriction that the pilot symbols 
must be chosen from the same signal set as the data symbols 
is removed, then the idea of data dependent pilot symbols 
can be extended easily to other modulation indices. It is also 
possible to use a longer pilot block, rather than a different set 
of channel symbols, to force the modulator to return peri- 
odically to the zero phase state. An illustrative example of 
the two alternatives is given later for the case of binary full 
response CPM with a modulation index of h=%. Also shown 
later is the proof that the proposed pilot symbol encoding 
technique can be generalized to all rational modulation 
indices of the form h=p/M. It should be pointed out that the 
"optimal** precoders for different integer p will have in 
general, different assignments of channel symbols to their 
transitions (but the basic trellis structure remains the same) 
Thus, a pilot symbol assisted M-level full response CPM 
schemes with h=l/M uses an even transmission frame size 
N with one data dependent pilot symbol per frame. The pilot 
symbol insertion rule is given by 



10 where 



(18) 



(19) 



15 



20 



25 



and the q^O's were defined in (2K4). Given mat q 0 (T>= 
q 2 (T) and that qoCr^iCD^^^l. the term 8 n 2 can be 
made independent of the pulse shape if we set c 1 =c 0 . In this 
case, 6 n _ 2 equal to nui/M modulo 2n when n is even, where 
m is an integer from the set {±1+3, . . . , ±(M-1)}. This 
means one can force 4>„ to zero by choosing c^^^-m. 
The above observation suggests that a pilot symbol assisted 
partial response CPM scheme with L=3 and h=l/M would 
use an even frame size N, with each frame consists of 
K-N-2 information symbols and 2 data dependent pilot 
symbols. The pilot symbol encoding rule is given by 



<**-i=**af=-*i if Qu<_2=+*rm/M mod 231, 



(20) 



30 



= ~m if 71 mod 2n 



(16) 



where c w is the pilot symbol of the k-th frame. Since this 
rule guarantees that <Jw=0 for any integer k. the output of the 
anti-aliasing filter at time kNT becomes 



35 



where c Jtv _ 1 and are the pilot symbols in the k-th frame. 
As in the case of full response signalling, the received 
sample at time kNT is given by (17) and is independent of 
the channel symbols. 

The general pilot symbol encoding rules for even and odd 
pulse lengths L are given in Table 1 below. 



TABLE I 



PUot Symbo l Encoding Rules. 



L N 
21 even or odd 



possible values 
of Q 



pilot symbols pilot encoding rules 



[kN-M "l 
c^, + 2 Z c. 



{±1, ±3, 



,±(2M-1)> 



mod4M 



21+ 1 



L n=-l+l J 



mod2M 



{±1,±3,...,*M-1)> c^. 



where, i = 0, 1 1-2 



where, i = 0, 1 1-1 



which is independent of any channel symbols. With appro- 
priate suppression of the noise terms n w/ by the channel 
estimator, these periodic pilot samples can provide accurate 
channel state information. 
Partial Response Signalling 

The idea of data dependent pilot symbols can be gener- 
alized to partial response CPM schemes with an arbitrary 
pulse length L. For illustration purpose, consider the case 
when L=3 and h=l/M. As in the case of full response 
signalling, the pilot symbol encoding technique reported in 



Following a siinilar convention used to define the phase 
states, the result of a modulo 2M (or a modulo 4M) operation 
in Table 1 can take on a negative value. As in the case of 
L=3, the rules listed in Table 1 are obtained by exploiting the 
60 symmetry in the pulse shape f(t). It was found that in 
general, the number of pilot synibols required per frame 
equal 

« JW-x-f 1 *", (21> 

I i+- L y Lodd 
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Compared to full response signalling, pilot symbol assisted 
partial response schemes require more pilot symbols. 
Fortunately, most of the CPM schemes of practical interest 
have a pulse length shorter than 6 symbol intervals. This 
means each frame will consist of at most 3 pilot symbols, 
which for a reasonable fade rate f d T. does not significantly 
affect the throughput of the modem. 

Nyquist-3 Frequency Modulation — N3FM 
For the class of partial response CPM schemes using a 
Nyquist-3 frequency pulse the term q^t) in (2) satisfies the 
criterion: 



10 



■I: 



L-l 
5 

0 otherwise 



(22) 



Although the received sample r 0 is independent of any 
data symbol, other elements in R are not Let the length of 
the frequency pulse f(t) be written in the form L=2k+j, where 
j=0 when L is even and j=l when L is odd. Also let (s^. . 
. . ,So^!, . . . .s^j) be the i-th legal waveform of length 2a 
symbols (including pilot symbols) that passes through the 
zero phase state at time t=0. Based on the pilot symbol 
encoding rules listed in Table L it can be shown that when 
a is less than L. then there are altogether J=M a+A such 
waveforms and they have the common characteristic mat 
s o <0 =A. When a is greater man L-L then the number of such 
waveforms will be increased by a factor of M 0 "^ 1 . Let the 
diagonal matrix S, be defined as: 



where it is assumed mat the pulse length L is an odd integer. 
Consequently, the term <t>„ in (5) becomes 



15 



St=dU2gis^aa, . . . t s^\ 
where i=0,L . . . J-l. Then it can be shown that 



(23) 



20 



where 6 = 



L-l 



Except for a time shift, the above equation is identical to the 
one derived earlier for full response CPM; see Equation 14. 
In other words, all N3FM schemes require only one pilot 
symbol per frame. 
The Channel Estimator and Demodulator 
The received samples in Equation 11 are passed to the 



25 



O f „ = (l/7)Vsr 
and 

y-i 

t=0 



(28) 



(29) 



(30) 



Note that is the covariance matrix of the column vector 

G=(g^E/- &*/)' complex gain terms, O^f^NJ is the 

covariance matrix of the noise vector n=( n -a/ no/)'* 1 *s 

an identity matrix of size 2oI+L and <b gC is the correlation 



channel estimator. The estimator consists of a front end filter 30 of go and G. The components of the and <t> gG matrices 



F and an interpolation filter H (see FIG. 3). The function of 
the front end filter is to suppress as much as possible the 
wideband noise in the received pilot samples in Equation 17. 
These filtered pilot samples are then interpolated to obtain 
fading gain estimates at other instants. 
The Optimal Front End Filter 

Without loss of generality, consider only the estimation of 
the fading gain at rime t=Q. i.e. g^ The estimated value is 
denoted by go and the estimation error is e=go-go- Note that 
the estimation error e is also complex Gaussian. 

In principle, one can estimate go directly from the 
received pilot sample r^ This results in a mean squared 
estimation error equal to the variance of the noise sample n^ 
As mentioned earlier, the noise variance equals f^N^ where 
f s -VT is the two-sided bandwidth of the anti-aliasing filter. 
In order to suppress the inters ymbol interference generated 
from the filtering operation, it is required that I be a large 
number. This translates into an unnecessary large mean 
squared error (MSE) if go is estimated directly from r^ 

To reduce the MSE in the fading gain estimates to below 
f,N 0 . a front end filter should be used to determine go from 
the block of received samples 

R=Kr^, ro.r, rjf, (25) 

where 2a is the span of the filter given in number of symbol 
intervals, and (*)t represents the transpose of a matrix. 
Specifically 

SofFR (26) 

From linear estimation theory, the optimal filter F, in the 
minimum MSE sense, is given by 



can be determined easily from the autocorrelation function 
of the fading process in Equation 9. Given the optimal front 
end filter in Equations 27-30. the mean squared estimation 
error is 



35 



(31) 



45 



where & Rg is the Hermitian transpose of & gR . The value of 
the MSE is used to design the interpolation filter. 

Consider the ideal case when the receiver has perfect 
knowledge about the pilot and data symbols at the frame 
boundary. In this particular case, the front end filter is given 
by 

where S 0 is the diagonal matrix constructed from the pilot 
and the data symbols at the frame boundary. The mean 
squared estimation error of this front end filter is 



50 



(33) 



55 



60 



which is independent of S 0 . Of course one expects that e 2 0 
is less man the original MSE e 2 in (31). To enjoy this 
potential reduction in MSE (and hence a reduction in bit 
error rate) in practice, one can adopt an iterative decoding 
approach (i.e. decision feedback, or multiple-pass) whereby 
S 0 is replaced by the corresponding decision from the 
previous decoding attempt. 
The Interpolation Filter 

The output of the front end filter at time kNT, k an integer, 
can be written as 



where 



wher e cfr g/? =(*/ 2) goR" is the correlation between go and R and 

0^ / H , /2)RR* is the covariance matrix of R. Note that the 65 

notation (•)* is used to represent the Hermitian transpose of 

a matrix. and 



(34) 



(35) 
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(36 > Numerical Results 

is the estimation error of the front end filter. The e A 's are iid w St™ 3 * 7 8 analytical and simulation results for 
zero mean complex Gaussian random variables with a MSK * noi malized fade rates, f 0 T, of 0.01 and 0 03 
variance of e 2 (31). The y^s are fed to the interpolation filter respectively, are shown. Also included in these figures are 
for the estimation of the fading process at other instants. For 5 P^onnanoe curves of MSK with perfect channel state 
reliable estimation, the frame size N must be chosen to , 0nmtlon (CSI > differential detection. In the 

satisfy the Nyquist sampling theorem. This means for a ^ the J*** fr ? m end mter <^rived previously is 

Doppler frequency of f^ 2Nf*T must be less than 1 In to K ^a ^^TTl^ W«>1, the frame size 
study, N is chosen to be approximately 1/(3^ tn^ ^^ f ° f V^ 03 frame si2e is N =^- n* 

loss of generality, one only coders the ^^SS i0 SS^*?^ ? 2 ^ ^ d *esizeofthe 

time zero and ends at time NT 31 ^polator is 10 (or p=5). The analytical results are 

vector Y. the esuma.orcompu.es observed at the output of our pilot symbol Le^^teX To 

the contrary, even though there is no interleaving, the MSK 

g .=H<«>y .„. ^ f^£ cxhib j? s 3 effect ~ m fast fadin g channels. 

fc <37) This can be seen fromFIG. 7 where the bit errorprobability 

where the optimal interpolation filter is given by decreases at a rate faster than the inverse signal-to-noise- 

ratio. ~' 



M *"" n • <»> As shown in FIG. 7. at a BER of 1(T 3 and a f D T=0 01 

Here «f„ is the covariance matrix of the gain vector X. e 2 I ?xJ 2 " P ? SS decoder Provides a 2.5 dB improvement in' 

is the covariance matrix of e. O^ 2 ! is the covariance „ , when compared to its one-pass counterpart, 

matrix of Y. and is the correlation of fc and X The (C) At " BER <* 10 a fad «= ™* of f o T=0. 0 3. the 

corresponding mean squared interpolation error is precoder 22 having the trellis structure shown in FIG 3 

^=o *-j> «•>» +e , rr'<t> ST"" 3 1 *^ B o COdin8 6aiD rclative to a memoryless 

' " »wt«*0 »«/<■) (39) mapper; see FIG. 8. 

where <t>(n)^ is the Hermitian transpose of <D(n) Note v> jj^ ^tJZJ 3 ^ wim BT=0 - 4 01 shown in 
that in the ideal case when the front end filter has side ? T an 1 d t 10 ^ The frequency pulse f(t) is truncated to L=3 
information about the data and pilot symbols at each frame ^ Ille * Tame size >«d interpolator size are identical 
boundary, the term e 2 in both (38) and (39) should be °„ li". ? MS ^ However - *e span of the premier is 
replaced by e 2 „ from (33). Finally, we like to point out that * 4 s >' mlx ' t ( Le - «=2>- The observations are 

the interpolation error u 2 „ is rather insensitive to the data „ * * ose MSK. In comparing the two sets of 

position n. 35 results - found that GMSK has almost identical perfor- 

The Vitcrbi Demodulator mance to MSK m terms of the BER. Although not shown 

The sequence of channel estimates {&,> are fed to a precoder 22 having the structure shown in FIG. 3 

coherent CPM demodulator implemented in the form of the P rov,des l a coding gain of 1.7 dB when compared to a 
Viterbi algorithm (VA). The VA selects the sampled wave- an mc i no ^ ,ess ma PP er - Similar results on GTFM can be found 

form M . . . ) i„ mc CPMttelhs toat J?*?' U ™ " there is presented the effect of pre- 

maximizes the conditional probabflity density function of v " 8 ^^JZ 1 Continuous Phase Frequency Shift 
the received samples {r„} given s and {&} Since me jf™.? (4CPFSK) system with a modulation index of h=W 
received samples are Gaussian and that the CPM signal has t * dccoder - The precoder 22 being used is the one 
a constant envelope, this is equivalent to selecting the „, "J^S *f trellis structure shown in FIG. 4. In the absence 
waveform whose metric OI Preceding, a Gray mapper was used instead. The frame 

size, the interpolator size, as well as the span of the front end 
M(')=r{r-^J. + ^.J.. f (40) niter are identical to those used in MSK. Also shown in FIG. 

U is toe performance curve of a differential detector with a 
is the largest. pre-detection filter given by Equation 27. As in the case of 

It was shown in above that the MSE (mean square errort m«?L- 7 )- * CTC * 00 imducible «ror floor and 

of the estimator can be significantly redded ffttto^L Provides about 1.7 dB improvement in SNR. 

infonnation about me pilofand data sSohl, toelra£ SmlS^f ^Sr****? ^^g Technique 
boundaries. This suggests an iterative <Cg be «n^^ encoding tochnique discussed above can 

which the decisions of the Viterbi demoduUtofa^edback ss rat,onal niodulation indices of the form 

to the channel estimator for the purpose of conscucune T P i "r™ and M is the °"mber of 

better fading gain estimates in the next derodff a^S ! leve ? s ; ^ additi °n. an illustrative example is 

Although some of the fed back il^Zf 2yT £ U al to ^Tp* 6 ^ ° f 3 moduJation ^ not 
erroneous, it is believed that as lone as the initial hit!™ ^ !! convemence in the discussion, only full 

rate (BER) is relatively low. toe impa^trf «ontou S X" «o T^f 8 

sion feedback is imnimal. As an example tte^LX • h= «f I i ad P « odd. then the phase state <t»„ for odd n 
mance of a 1-pass (ie. no decision^dbaST) Triors" r JZZl ^ ^ l ° force * e to 

and a 3-pass decoder wfll be compared P T?T ^ ^? phase stote 81 ^ ( D+1 ) T - we must be able 

Despite the occasional erroneous decision feedback, the toarlo^^L ^ T ' ' ' such 

interpolator used in the second pass (and un> ofTmniTir^ « that cpjt/M=-mn/M modulo 2jc. This is equivalent to finding 
pa^coder is obtained^SS^ ^l^sfby ^3 '^"t^.^T 0 m0dul0 2M ' where c^ 
the term e 2 0 from (33). Wby m-m+M. and both c and m are odd integers from the set 

iU.l. . . . 2M-1}. It is easy to show that such a solution 
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always exists. Consequently, one can always find a valid continuous phase modulation scheme and means for decod- 
symbol from the signal alphabet to force the modulator to ing the demodulated signal to recover the estimated trans- 
return to the zero phase at even symbol time. A similar mitted data. 

analysis can be used to show that for an even integer p. it is 3. Apparatus according to claim 1, wherein said means for 

also possible to find a valid channel symbol to return the 5 receiving and filtering includes an anti-aliasing filter 
modulator to the zero phase state in one step. Thus, it can be designed to miiumize the total intersymbol and adjacent 
concluded that the proposed pilot symbol encoding tech- channel interferences. 

nique can be used in all M-level CPM systems with a 4. Apparatus according to claim 2. wherein the trellis 

rational modulation index of h=rVM. structure is that of a full response continuous phase modu- 

For modulation indices not equal to p/M, the pilot symbol 10 larion scheme within an identical number of modulation 
encoding rules must be tailored to the particular value under levels and modulation index to that of the underlying 
consideration. As an illustrative example, we consider the continuous phase modulation scheme, 
case of a full response binary CPM schemes with a modu- 5. Apparatus according to claim 2, wherein transitions of 

lation index of h=%. There are two possible sets of pilot the modulator joining at the same state are assigned different 
symbol encoding rules and they are discussed as follows. 15 channel symbols, transitions emanating from the same state 

When h=%. then the phase state at time nT. n odd, is of are assigned different channel symbols, transitions joining at 
the form the samG m * e fell** nave identical information 

content and transitions emanating from the same state in the 
mn (51) trellis have their information content gray-coded in accor- 

A 20 dance with the squared Euclidean distances between wave- 

where m e {±143}. Thus we can force ^ to zero by form pairs associatedwith that state 
making c^=Ui a data dependent pilot symbol: In this case, 6. Ar*aratus according to claim 4 ^herein 

one is in~ sence using a 4-level modulation for the pilot joining at the same state are assigned different channel 
syrnbols. Some degradation in error performance is exo^ symbols, transitions emoting from the same state arc 
ulL 25 assigned different channel symbols, transitions joining at the 

^Alternatively, one can use three pilot symbols c n+1 . Cjt+2 . same state in the trellis have identical Morrnation content. 
c n+ 7 Return me modulator to the Vero phaTstate^at time transitions emanating from the « ^^^^ 
(n+W. The pilot symbol encoding rules in this case are their information content gray-coded according to the 
1 c =c J* - 1 if m=l squared Euclidean distances between waveform pairs asso- 

7 r ""' f ~ 2 II n+r l7 ^ ' 30 dated with that state. 

T ^Q^™Ciifm=3 7. Apparatus according to claim 1, wherein said pilot 

* ^'^ 2 =c n+3 i tf ^-3' symbol insertion device uses at least one pilot symbol for 

This%proih Je^uires 3 pilot' symbols per frame and hence every N transmitted symbols to force the modulator to return 
leadsTo a drop in the throughput of the modem However. periodically to a set of known phase states, 
both the data°L pilot symbSl/are chosen from the same set 35 « . ^tus according \^ 7 ^^*^ U ?f 
of {±1} symbols the pilot symbols may be substantially different from those 

Accordingly, while this invention has been described with used for me data symbols. 
rcfercnceTmustrative embodiments, this description is not » Apparatus according to claim 1. wherein said channe 

intended to be construed in a limiting sense. Various modi- estimator provides periodic, minimum mean channel 
fixations of the illustrative embodiments, as well as other <o estimates and being designed ™*™^£™*°2Z 
embodiments of the invention, will be apparent to persons more of *> u ^ 
skilled in the art upon reference to this description. It is to noise ratio, the Doppler frequency, die 
therefore contemplated that the appended claims will cover function of the channel fading process. f 6 "^ 00 
any such modulations or embodiments as fall within the of the pilot symbols and the set of P 055 ^^^^^ 
true scope of the invention. of a g* vcn duratlon mat ^ ata me pilot symbo1 

We claim- insertion point. 

1. Apparatus for transmitting and receiving digital infor- 10. Ar^ararus acccctog to claim 1. wher ein ^dd^nd 
mation over mobile communication channels, said apparatus estimator mcludes a front end filter and an int erpo aUon 

... filter, said front end filter providing periodic, minimum 

comprising. . . . 50 mean square channel estimates and being designed in accor- 

(a) a pilot symbol insertion device for periomcally insert- * more of the foUowing channel and system 
ing data dependent pilot symbols into frames of digital oonditfonK the signal to noise ^c, Doppler frequency, 
data for subsequent channel estimation; ^ artoconeUtion function of the channel fading process. 

(b) a continuous phase modulator for modulating and me sym bol insertion points and the set of possible 
transmitting over a mobile communication channel the 55 ^j^uous phase modulation waveforms of a given duration 
frames of digital data, wherein said pilot symbol inser- thajt ccnter around the pilot symbol insertion point 

tion device uses data dependent pilot symbols to force ^ Apparatus according to claim 10. wherein the inter- 

the continuous phase modulator to return periodically portion filter when coupled with said front end filter pro- 
to a set of known phase states; vi<tes channel estimates at all desired time instants. 

(c) means for receiving and filtering the transmitted data; ^ 12. Apparatus according to claim 11, wherein said inter- 

(d) a channel estimator for estimating the channel ampli- polation filter is a digital filter the design of which is based 
tude and phase distortion from the received pilot signal upon the signal to estimation error ratio at the output of the 
at different time instants; and front end filter, the Doppler frequency and the autocorrela- 

(e) means for demodulating the received signal. tion function of the fading process. 

2. Apparatus according to claim 1. including a precoder 65 13. Apparatus according to claim 1. wherein said demodu- 
employing a finite state preceding method whose operation lator means is a coherent continuous phase modulation 
is based upon a trellis structure substantially the same as the demodulator implemented in the form of a Viterbi algorithm 
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which operates on the received signal and the channel r« a ~. . • 

estimates provided by the channel eslLatoT Jl^T* Si * nal bascd "J™ 

14. Apparatusaccordngtoclaim 1. wherein said demodu- anT * ^ estimator ". 

^£^7^°^ (8) decoding the demodulated signal 

-5- Apparatus according to claim 1. wherein said demodu- ^T^T^^TSZ n^hf 

latormeansisaniultiple-passjoimdatadete^andchannel ,aUon scbcme ~ continuous phase modu- 

estimator by which temporary decisions made by a data ln 2i A method according to claim 20 wherein said trellis 

detector in the current demodulation pass arc used to refine f^Tl ** ° f 8 ltal res P onse «»""™ous phase modu- 

the channel estimates in the next demodulation pass S< ^V°^ ™?> identical number of modulation levels 

16. Apparatus according to claim IS. wherein said ^ modu r a * oa lndex t0 «>at of the underlying continuous 
demodulator means utilizes a Viterbi algorimm ^f^" HZ* 

17. Apparatus according to claim IS wherein said « of ^ n^T.^f * c F < ! r * n 8 to claim 21. wherein transitions 

.cans utilizes a reduced comple^^ St"32T » 15= emTat^g £f£ 

IS. Apparatus according to claim 15. wherein said mul- X^JSZ2£l ^2S£S 

uple channel pin, data detector and channel estimator M formation content and transitions enmatog from £ 

uidudes a front end channel estimation filter mat provides " same state * *e trellis have their infonnafon intent 

periodic, minimum mean squared channel estimates bascd eroded in accordance with the squared Euclidean dis 

ZZ^J? , 10 DOiSe ^ D °W ,ler fr«!»«cy. the ^ between wavefonn pairs associated with that state, 

autocorrelation function of the channel fading process, the ™' A m f^° d acco ««"g 'o claim 20. wherein at least one 

mseruon points of the pilot symbols and me continuous 25 T*?*?. USed for N transmitted symbols to 

phase modulation waveform of a given duration that center £ f modulalor "> return periodically to a set of known 
around the pilot symbol insertion point and specified by the 

previous demodulation attempt 25. A method according to claim 20. wherein said channel 

19. Apparatus according to claim 15. including an inter "bating step provides periodic, minimum mean square 

polation filter that uses the channel estimates provided bv 30 7 estimates in accordance with one or more of the 

the front end channel estimation filter T^ZZ^l S TS. ^ me ««"»- 

(a) encoding a stream of digital data including inserting ^ . **' A mctbod according to claim 20. wherein said 
periodically data dependent pilot symbols into frames ^T^ 1 " 8 Step utmzes a ^Uerbi algorithm that operates 
of the digital data; on the received signal and the channel estimates 

(b) performing continuous phase modulation on the demodl^fl^l according to claim 20. wherein said 
frames of data, wherein inserting the data dependent 2Z£^£f. ^ * redUCed ««»P»e«ty sequence 
pilot symbols forces the continues phj fJSSSl « San^ esSnateT" °° Md * e 
to return periodically to a set of known phase states- 20 A nwK,J" 

(d) detecting the transmitted data; ^ b X ^ data detector in the current demodulation pass 

(e) estimating from the received pilot signal the channel's ^ * refine 1113,11161 estimates ^ the next demodu- 
amplitude and phase distortion at different time ***** 



instants; 
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